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RESUMEN 

El resumen sera traducido al espanol por los editores. We study the formation and evolution of self- 
interacting dark matter (SIDM) halos. We find analytical, fully cosmological similarity solutions for their 
dynamics, which take proper account of the collisional interaction of SIDM particles, based on a fluid approx- 
imation derived from the Boltzmann equation. These similarity solutions are relevant to galactic and cluster 
halo formation in the cold dark matter (CDM) model. Different solutions arise for different values of the 
dimcnsionlcss collisionality parameter, Q = ap\,r a . where a is the SIDM particle scattering cross section per 
unit mass, pb is the cosmic mean density, and r s is the shock radius. For all solutions, a flat-density, isothermal 
core is present which grows in size as a fixed fraction of r s , pumped by cosmological infall. Accordingly, core 
collapse must in general be delayed until infall becomes negligible, contrary to previous analyses based on iso- 
lated halos, which predict core collapse in a Hubble time. Our solutions agree with N-body simulations, which 
match observed galactic rotation curves if a ~ [0.56 — 5.6] cm 2 g _1 , implying Q ~ [6.2 x 10~ 7 — 3.6 x 10~ 5 ]. 



Similar profiles also arise for Q 
not previously simulated. 



[1.37 x 10~ 2 - 1.7 x 10" 1 ], or a ~ [1.2 x 10 4 - 2.7 x 10 4 ] cm 2 ff -\ a regime 



ABSTRACT 



We study the formation and evolution of self-interacting dark matter (SIDM) halos. We find analytical, fully 
cosmological similarity solutions for their dynamics, which take proper account of the collisional interaction 
of SIDM particles, based on a fluid approximation derived from the Boltzmann equation. These similarity 
solutions are relevant to galactic and cluster halo formation in the cold dark matter (CDM) model. Different 
solutions arise for different values of the dimcnsionlcss collisionality parameter, Q = api,r s , where a is the 
SIDM particle scattering cross section per unit mass, pb is the cosmic mean density, and r s is the shock 
radius. For all solutions, a flat-density, isothermal core is present which grows in size as a fixed fraction of 
r s , pumped by cosmological infall. Accordingly, core collapse must in general be delayed until infall becomes 
negligible, contrary to previous analyses based on isolated halos, which predict core collapse in a Hubble time. 
Our solutions agree with N-body simulations, which match observed galactic rotation curves if a ~ [0.56 — 
5.6] cto 2 ^- 1 , implying Q ~ [6.2 x 10" 7 -3.6x 10~ 5 ]. Similar profiles also arise for Q ~ [1.37 x 10" 2 -1.7x 10" 1 ], 
or a ~ [1.2 x 10 4 — 2.7 x 10 4 ] cm 2 g~ l , a regime not previously simulated. 

Key Words: COSMOLOGY:DARK MATTER — COSMOLOGY: LARGE-SCALE STRUCTURE OF 
UNIVERSE — GALAXIES: KINEMATICS AND DYNAMICS 



1. INTRODUCTION 

The self-interacting dark matter (SIDM) model 
was proposed by Spergel & Steinhardt (2000) to re- 
solve the discrepancy between the inner density pro- 
file of cold dark matter (CDM) halos found by N- 
body simulations and observations of galactic rota- 
tion curves. N-body simulations have a density cusp 
(p oc r 13 where (3 ranges from -1 (Navarro, Frenk 
& White 1997; hereafter "NFW" ) to -1.5 (Moore 
et al. 1999; hereafter "Moore profile")), while ob- 
servations of dark-matter dominated dwarf and LSB 
galaxies indicate flat-density (soft) cores. SIDM par- 
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tides undergo nongravitational, microscopic interac- 
tion which is strong enough to produce a soft core by 
heat conduction. Cosmological N-body simulations 
which incorporate a finite scattering cross-section a 
for the SIDM particles show that this scheme suc- 
cessfully produces soft cores in halos (e.g. Dave et 
al. 2001; Yoshida et al. 2000b). 

We study the formation and evolution of SIDM 
halos analytically with a proper treatment of cosmo- 
logical infall (see also Ahn & Shapiro 2002). Pre- 
vious analytical studies (Balberg, Shapiro & Ina- 
gaki 2002) and N-body simulations of isolated halos 
(Burkert 2000; Kochanek & White 2000) neglected 
cosmological infall which might delay the core col- 
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lapse. We find that, for a realistic range of infall 
rates, the collapse of SIDM cores can be completely 
inhibited. Our results agree with N-body simulations 
for nonisolated halos, but also find that a range of 
higher scattering cross sections than previously sim- 
ulated can also produce soft cores. 

2. SELF-SIMILAR EVOLUTION OF SIDM 
HALOS 

We use a fluid approximation derived from the 
Boltzmann equation to describe the dynamics of 
SIDM halos. A spherically symmetric, arbitrarily 
collisional system with isotropic random motions can 
be described by fluid equations for an ideal gas with 
7 = 5/3 (Bettwieser 1983). 

For a collisionless system, collapse is "adiabatic" 
(i.e. no conduction), and if the initial overdensity 
has a scale-free power-law form, 

where e is a positive constant, the object which grows 
by adiabatic infall is self-similar (Fillmore & Gol- 
dreich 1984). Heat conduction generally breaks self- 
similarity by introducing an additional length (time) 
scale, but when r s oc t 2 , or e = 1/6, self-similarity 
is preserved even in the presence of heat conduction 
by SIDM particles with constant a. 

3. SELF-SIMILAR CDM HALOS (e = 1/6; NO 
CONDUCTION) 

The condition which is required to preserve self- 
similarity of the dynamics of SIDM halos deserves 
separate attention. We find that the adiabatic solu- 
tion with e = 1/6 is well fit, for example, by an NFW 
profile with concentration parameter 3 < cnfw _ 4. 
The inner density slope is —1.27, which is between 
-1 (NFW) and -1.5 (Moore profile) (see Fig. 1). 

The theory of structure formation from density 
peaks in the Gaussian random noise distribution of 
initial density fluctuations gives an interesting clue 
to this correspondence. According to Hoffman & 
Shaham (1985), for a power-law power spectrum 
of initial fluctuations, P(k) cx k n , the initial den- 
sity profile of density peaks can be approximated 
as oc r K , where k = 3e = n + 3. The value 
e = 1/6 corresponds to re = 1/2, or n = —2.5: 
n = —2.5 roughly corresponds to galaxy-mass struc- 
tures in ACDM. 

4. SELF-SIMILAR SIDM HALOS (e = 1/6; WITH 
CONDUCTION) 

We find that different similarity solutions of 
SIDM halos arise for different values of the collision- 
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Fig. 1. Halo density profile for adiabatic solution with 
e = 1/6 for standard CDM halos, compared to NFW 
and Moore profiles. The adiabatic solution, whose inner 
slope is -1.27, is a good fit to density profiles of CDM 
halos from N-body simulations. 

ality parameter, Q = apt,r s . High Q means high col- 
lision rate. We also find that there are two different 
regimes: the low-Q regime (Q < Qth = 7.35 x 1CU 4 ) 
and the high-Q regime (Q > Qth)- 

In the low-Q regime, the collision mean free path 
in the core is larger than the size of the halo. In 
this regime, the flattening of the core density profile 
increases as Q increases (see Fig. 2). 

In the high-Q regime, the mean free path in the 
core region is smaller than the size of the halo. The 
flattening of the core density profile decreases as 
Q increases, because the mean free path decreases 
(diffusion limit; see Fig. 2). The limiting case of 
Q = oo corresponds to the solution with no con- 
duction, which agrees with N-body simulations with 
infinite cross-section (Yoshida et al. 2000a; Moore 
et al. 2000). In this case, the nonadiabatic solution 
approaches the adiabatic solution. 

Our similarity solutions naturally explain how 
cosmological infall affects SIDM halos dynamically. 
New matter which falls into the halo from outside 
gets shock-hcatcd (i.e. virializcd) to such an extent 
that the temperature is always greater in the halo 
than in the core. Accordingly, in our similarity solu- 
tions, core collapse is completely prohibited. There- 
fore, core collapse will be inhibited in general until 
the infall rate becomes negligible at late times in the 
evolution of individual halos, as reported for N-body 
simulations by Wechsler et al. (2002). 

5. THE ALLOWED RANGE OF SCATTERING 
CROSS-SECTION FOR AN SIDM UNIVERSE 

We find that a relatively narrow range of Q val- 
ues, 

[0.62-3.6] x 10- (A;) (jj^r), (2) 
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Fig. 2. Similarity solution dimensionless profiles for low-Q (left panels) and high-Q (right panels) regimes. Even though 
values of Q are very different, profiles of high-Q solutions are indistinguishable from those of low-Q solutions. 



characterizes the entire range of halo masses from 
10 6 /i _1 M Q to lO 13 /!" 1 !^, in the currently-favored 
ACDM universe (see Fig. 3). 

The range a = [0.56 — 5.6] cm 2 g _1 is the pre- 
ferred range of the scattering cross section found by 
cosmological N-body simulations for the ACDM uni- 
verse to match observed galactic rotation curves (e.g. 
Dave et al. 2001). Equation © with h = 0.7 then 
yields Q = [6.2 x 10" 7 - 3.6 x 10" 5 ], which is in the 
low-Q regime. 

As shown in § 0J there also exist high-Q so- 
lutions which yield profiles which are quite simi- 
lar to the low-Q solutions which produce observa- 
tionally acceptable soft cores. We find that Q = 
[1.37 x 10~ 2 - 1.7 x 10~ 4 ] in the high-Q regime pro- 
duces profiles similar to those in the low-Q regime for 
Q = [6.2 x 10~ 7 - 3.6 x 10" 5 ]. From the relation- 




Fig. 3. The collisionality parameter Q vs. mass of halos 
at their typical formation epoch for a = 5.6cm 2 g -1 and 
h = 0.7 in ACDM universe. 



ship between a and Q (equation |2J), we predict that 
a = [1.2xl0 4 — 2.7xl0 4 ] cm 2 g~ 1 can also produce ac- 
ceptable soft cores, therefore. Cosmological N-body 
simulations which incorporate SIDM particles with 
a ~ [5 x 10 3 — 5 x 10 4 ] cto 2 5 _1 would, therefore, be 
of interest. 

This work was partially supported by NASA 
ATP Grant NAG5-10825 and Texas Advanced Re- 
search Program Grant 3658-0624-1999. 
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